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Abstract 

Recently, a holographic argument was used to relate the saturation 
exponent, jbfkl, of four-dimensional Yang-Mills theory in the Regge 
limit to the Choptuik critical scaling exponent, 75^, in 5-dimensional 
black hole formation via scalar field collapse pp. Remarkably, the nu- 
merical value of the former agreed quite well with previous calculations 
of the latter. We present new results of an improved calculation of 
75d with substantially decreased numerical error. Our current result 
is 75d = 0.4131 ± 0.0001, which is close to, but not in strict agreement 
with, the value of jbfkl = 0.409552 quoted in pQ. 

The critical behaviour exhibited by gravitational collapse to black holes, 
first discovered by Choptuik [2] (see pHH] for reviews), is a small scale, strong 
gravitational field effect that is by now well understood at the qualitative 
level. As well, key physical quantities such as the associated critical solution 
and scaling exponents have been calculated numerically. Precise analytic 
derivations of these quantities are still, to the best of our knowledge, lacking. 
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Recently, intriguing holographic arguments were used to relate the strong 
field gravitational scaling behaviour in spherically symmetric scalar field col- 
lapse in five space-time dimensions to scaling in a suitable, 4-D weakly cou- 
pled gauge theory pQ. Specifically, it was conjectured that exponential growth 
of the scattering amplitude for two hadrons in the one-pomeron exchange 
approximation corresponded via a holographic mapping to the exponential 
growth away from the critical solution that leads to Choptuik scaling in five 
dimensions. Thus, it was argued, the BFKL scaling exponent Jbfkl should 
equal the 5-dimensional Choptuik scaling exponent 75^. 

As shown by the authors of PQ, a numerical calculation of the BFKL 
scaling exponent yields the value Jbfkl = 0.409552. On the other hand, the 
most accurate previous calculation yields a 5-D Choptuik scaling exponent of: 
75d = 0.412±0.004 [5J. A somewhat lower result , namely 0.408, was reported 
in [6], but with an error of ±0.008. The apparent agreement between •Jbfkl 
and provides startling support for the conjecture relating 5d gravity to 
pomeron exchange in 4d Yang-Mills theory, and for the general validity of 
the holography hypothesis. 

The purpose of this note is to present the results of a more accurate cal- 
culation of 75^ and the related echoing period A 5 ^. Recall that the Choptuik 
scaling exponent is most directly observed in the simple scaling law obeyed 
by the horizon radius near criticality [2] 

Rah oc (a - a*) 7 . (1) 

where Rah is the horizon radius on formation, a is a parameter describing 
the initial data whose critical value a* separates initial data that lead to 
blackhole formation from those that lead to dispersal of the scalar field. The 
echoing period A 5 d derives from the discrete self-similarity of the critical 
solution, and can be extracted from the periodicity of the critical scalar field 
solution: 

XM = X (re A ™,te A ™) (2) 
The new values we obtain for these quantities are 

7 M = 0.4131 ±0.0001 A 5d = 3.225 ± 0.003. (3) 

These results are consistent within error with previous calculations, but 75^ 
no longer agrees precisely with ^bfkl- Nevertheless, given the uncertainties 
in the calculation of 'Jbfkl, a discrepancy of this order (less than 1%) is 
perhaps not surprising [7]. 
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We now describe how was obtained. Our formalism and code are 
basically the same as that used in [5], which in turn was based on the analysis 
of [8]. A couple of crucial changes have permitted the increase in the accuracy 
by a factor of about twenty. 

The calculation can generically be formulated in terms of Einstein gravity 
in D space-time dimensions minimally coupled to a massless scalar field [8]: 



- j d D x ^m\d x ?. (4) 

Spherical symmetry is imposed by requiring: 

ds 2 D \ = g a pdx a dx 13 + r 2 (x a )dVL( yD _ 2 ), 

X = X(x a ) (5) 

where dQ^o-2) is the metric on S D ~ 2 and a, (3 = 1,2. A relatively simple 
set of equations can be obtained by the following field redefinition, which is 
motivated by 2-dimensional dilaton gravity [8j [9] : 

D-2 fr\ D ~ 2 

(6) 



r ' 8(D-3) \l 
9*p := 4> {D ~ mD - 2) 9«p, (7) 

where r is the optical scalar and (f) is proportional to the area of a {D — 2)- 
sphere at fixed radius r. 

We go to double null coordinates first introduced by Garfinkle [TO] : 



ds 2 = —2lg(u,v)(j)'(u,v)dudv (8) 
in which the relevant field equations take the form: 



D — 6 / r \ d-2 



(9) 



^ = 32vr( X ') 2 (10) 

(<t>x'y + m' = o, (ii) 
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where the prime and dot denote v and u derivatives, respectively. We hence- 
forth restrict attention to five dimensions by setting D = 5, in which case 



a/3 



The general method we use starts with the choice of a one parameter 
family of initial data for the scalar field x- Coordinate invariance allows 
the initial value of <ft to be chosen arbitrarily, while the metric component g 
is obtained by integrating the constraint Eq. ffTUj) . We then evolve the data 
using the remaining dynamical equations (Q and (TTTj) until either a horizon 
forms, or the field disperses. Horizon formation is signalled by the vanishing 
of the null expansion scalar, which for spherically symmetric geometries is 
given by: 

AH = gffidMt. (12) 



This calculation is repeated for different values of some parameter, a, describ- 
ing the initial data in order to determine the critical value, a*, separating 
black hole final states from dispersed final states. Once a* is determined, we 
calculate the horizon radius at formation, Rah, for a wide range of super- 
critical values of a near a* and then plot \d.(Rah) vs ln(a — a*). This gener- 
ically yields an (almost) straight line (see Figure [T]) from which the scaling 
exponent and echoing period can be obtained by fitting the ln-ln plot to the 
following 5 parameter function: 



In (R A h) = c + 7sd ln(a - a*) + c x sin 



27r \, / \ 
— in(a - a*) + c 2 

05d/ 



(13) 



where {co,ci,C2} are family dependent constants. The periodic deviations 
from a line in the ln-ln plot yield the echoing period via the relation A§d = 
^l5d^5d [llj- A more accurate value for A can be obtained by measuring 
the period of oscillations of the matter field at the origin as a function of 
—ln(u* — u), where u* is the a = a* value of the coordinate u at collapse (See 
Figure [2]). 

Our specific numerical scheme used a v ('space') discretization to obtain 
a set of coupled ODEs: 



h(u, v) — > hi(u) 



(U,V) 



(j)Au) 



(14) 



where i — 0, . . . , (N— 1) specifies the v grid and h = x+20x'/0'. Initial data 
for these two functions were prescribed on a constant v slice, from which the 
function g(u,v) was constructed. 
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Figure 1: Plot of Iti(Rah) vs ln(a — a*) (top graph). The value of a* was 
adjusted within a small range in order to increase the quality of the fit and to 
determine the error in 75^. In the bottom graph, the linear portion of the fit 
is subtracted from Iti(Rah)- The residuals are well described by a sine wave 
(solid curve) and clearly show the oscillation expected due to the discrete 
self-similarity of the critical solution. 
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Figure 2: Plot of the matter field at the origin vs —ln(u* — u). The value of 
was adjusted within a small range to increase the quality of the fit and to 
determine the uncertainty of the best estimate of the self-similarity period. 
A 5 d is determined by calculating the period of the oscillation (the solid line 
represents the best fit to the data.) 



6 



In all cases we started with Gaussian initial data 



x(« = O,0) = a0 2/3 



cxp 



a 



(15) 



The initial scalar field configuration \{ u — 0, <p) is most conveniently specified 
as a function of rather than r. (Recall that oc r 3 .) The value of the dilaton 
field at each grid point on the initial surface was chosen for simplicity to be: 



H(u = 0,v) = (i + l)Av, 



(16) 



where Ad is the initial grid spacing. The initial values of the other functions 
were determined in terms of the above by computing the integrals gi using 
Euler's method for equally spaced points. 

The results in Eq.(|3]) were obtained with an initial spacing Av = 0.0012 
along with an initial grid size of N = 2048. We further verified these results 
by performing calculations at several smaller grid sizes (with proportionately 
larger initial grid spacings) (See Tabled]). To the stated error, the scaling 
exponent was independent of initial grid size, while echoing period A^d varied 
between about 3.223 and 3.226. This determined the error bars on As,i given 
in Eq.©. 



Initial Grid Size 


Initial Grid Spacing 


75d 


A 5d 


512 


0.0048 


0.4130 ±0.0001 


3.2264 ±0.0005 


800 


0.003072 


0.4130 ±0.0001 


3.2230 ±0.0006 


1024 


0.0024 


0.4131 ±0.0001 


3.2223 ±0.0006 


1536 


0.0016 


0.4131 ±0.0001 


3.2264 ±0.0004 


2048 


0.0012 


0.4131 ±0.0001 


3.2246 ±0.0005 



Table 1: Choptuik 5-d critical exponent and echoing period calculated with 
varying initial grid sizes. 

Evolution in the 'time' variable u was performed using the 4 th order 
Runge-Kutta method. The general scheme is similar to that used in [12], 
together with some refinements used in [TO]- This procedure was also used 
for the 3— dimensional collapse calculations in |13j . The accuracy of the 
calculation was improved by decreasing the time steps as the calculation 
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progressed according to the formula: 

Am = Auq ( 1 - 0.999 exp - 



( 



w* — u 



0.1 



) 



2 



) 



(17) 



where Au = 0.001 is the initial time step spacing and is the estimated 
value of the coordinate u of the critical solution. Further improvements in 
accuracy were made by decreasing the value of Aw periodically during the 
calculation, or once apparent horizon formation commenced, i.e., once the 
function AH (</>) developed an extremum. Time steps ranged from Am = 10 -3 
at the beginning of the calculation to about Au = 10 -11 near the end. 
The boundary conditions at fixed u are 



where k is the index corresponding to the position of the origin = 0. All grid 
points < i < k — 1 correspond to ingoing rays that have reached the origin 
and are dropped from the grid. The boundary conditions are equivalent to 
r(u,u) = 0, g\ r =o = g{u,u) = 1. Notice that for our initial data, and 
hence are initially zero, and therefore remain zero at the origin because 



The most significant improvement in the numerical code over [5] , involved 
systematically decreasing the grid spacing Av using a method similar to that 
of [10]. At regular intervals in ln(u* — u), the exterior half of the gridpoints 
were discarded. The number of gridpoints was then restored by using a 
4 th order Lagrange interpolation between the remaining grid points, thereby 
doubling the resolution. Typically, by the end of the calculation the grid 
spacing was decreased to approximately 10 -6 . This allowed the accuracy of 
the horizon position to be maintained over the entire ln-ln plot and reduced 
significantly the overall uncertainty in 75^. 

At each u step, a check was made to see if an apparent horizon had formed 
by observing the function AH in Eq. (|T2|) whose vanishing signals the forma- 
tion of an apparent horizon. For each run of the code with fixed parameters 
a, a, and <po, this function was scanned from larger to smaller radial values 
after each Runge-Kutta iteration for the presence of an apparent horizon. In 
the sub-critical case, all the radial grid points reached zero without detection 
of an apparent horizon, which signals pulse reflection. In the super-critical 
case, black hole formation was signalled by the vanishing of AH at some 
finite radius. Since the code crashes when AH < 0, so the the time evo- 
lution was slowed down as AH — > + in order to obtain the most accurate 



0fc = g k = 1, 



(18) 



of Eq.jnj. 
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estimate of the apparent horizon position. For any given run, the apparent 
horizon position (Rmin) was assigned the value of the grid point for which 
AH reached a minimum during the last surviving iteration. 

For our initial runs we fixed the location and width of the Gaussian initial 
data to be 0o = 1 and a = 0.3, respectively. We then varied the amplitude 
a to determine the critical amplitude a* via a binary search. The code was 
written using double precision in C thus, the binary search was terminated 
when a* was determined to 15 digits. We then did runs in the super-critical 
region a > a* to determine the dependence of the horizon radius Rah on the 
amplitude. It was found that numerical error would introduce noise in the 
ln-ln plot for the runs closest to a* and so only super-critical runs for which 
(a — a*) /a* > 1CT 13 were fitted to determine 75^. To improve the quality 
of the fit Eq. (fT3l) . the value of a* was adjusted within the range of values 
between super-critical and sub-critical collapse obtained during the binary 
search. The error associated with 75^ represents the range of fit results when 
a* is adjusted within that range. The universality of the critical solution 
was verified by starting runs a = a* and varying the width and center of the 
initial data profile. As shown in Table [U the results for 75^ agreed within 
error in all cases. 



Critical parameter 


75d 


a 


0.4131 ±0.0001 


a 


0.4131 ±0.0003 


00 


0.4133 ±0.0004 



Table 2: Choptuik scaling exponent in 5 dimensional super-critical collapse. 
We verified the universality of the critical solution by varying the width and 
center of the initial data profile. The scaling constants calculated by varying 
a and 0o are in agreement with the results obtained by varying the amplitude 
of the initial pulse. Universality was tested at an initial grid size of 2048. 

To summarize, we have obtained fairly accurate values for the Choptuik 
scaling exponent and echoing period in spherically symmetric scalar field 
collapse in five dimensions: 

7 5rf = 0.4131 ± 0.0001 A 5d = 3.225 ±0.003. (19) 

These results agree within error with all previous calculations, but the value 
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of 75^ does differ somewhat from 'Jbfkl- The discrepancy is small but po- 
tentially significant. It should, however, be kept in mind that the value 
calculated in [T] is strictly speaking exact only in the asymptotic high energy 
region and it is not yet completely clear where this region begins. Hence there 
may be corrections due to the usual running of the strong coupling constant, 
saturation and other effects [7]. Thus although the proposed holographic 
connection between Jbfkl and 75^ is compelling, more work is required in 
order to confirm its validity. 
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